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ABSTRACT 

We present high precision radial velocities (RVs) of double-lined spectroscopic binary stars HD78418, 
HD123999, HD160922, HD200077 and HD210027. They were obtained based on the high resolu- 
tion echelle spectra collected with the Keck I/Hires, Shane/CAT/Hamspec and TNG/Sarge tele- 
scopes/spectrographs over the years 2003-2008 as a part of TATOOINE search for circumbinary 
planets. The RVs were computed using our novel iodine cell technique for double-line binary stars 
which relies on tomographically disentangled spectra of the components of the binaries. The precision 
of the RVs is of the order of 1-10 m s~^ and to properly model such measurements one needs to 
account for the light time effect within the binary's orbit, the relativistic effects and the RV variations 
due to the tidal distortions of the components of the binaries. With such proper modeling, our RVs 
combined with the archival visibility measurements from the Palomar Testbed Interferometer allow 
us to derive very precise spectroscopic/astrometric orbital and physical parameters of the binaries. 
In particular, we derive the masses, the absolute K and H band magnitudes and the parallaxes. The 
masses together with the absolute magnitudes in the K and H bands enable us to estimate the ages 
of the binaries. 

These RVs allow us to obtain some of the most accurate mass determinations of binary stars. The 
fractional accuracy in msini only and hence based on the RVs alone ranges from 0.02% to 0.42%. 
When combined with the PTI astrometry, the fractional accuracy in the masses ranges in the three 
best cases from 0.06% to 0.5%. Among them, the masses of HD210027 components rival in precision 
the mass determination of the components of the relativistic double pulsar system PSR J0737-3039. 
In the near future, for double-lined eclipsing binary stars we expect to derive masses with a fractional 
accuracy of the order of up to ^0.001% with our technique. This level of precision is an order of 
magnitude higher than of the most accurate mass determination for a body outside the Solar System 
— the double neutron star system PSR B1913+16. 

Subject headings: binaries: spectroscopic — stars: fundamental parameters — stars: individual 
(HD78418, HD123999, HD160922, HD200077, HD210027) — techniques: radial 
velocities 



1. INTRODUCTION 

The first observations of a spectroscopic binary star, C UMa (Mizar), were announced by Edward C. Pickering 
(1846 -1919) on 13 November 1889 during a meeting of the National Academy of Sciences in Philadelphia (jPickeringl 
|1890| ). A similar announcement about /3 Per (Algol) was made by Herm an C. Vogel ( 1841-1907) on 28 November 1889 
during a session of Konglich-Preussiche Akademie der Wissenschaften (Vogel 1890a, b). Even though a transatlantic 
telegraph cable had been available since the 1860s, it was quite unusual timing for a pre-astro-ph era. Pickering noted 
that the K line of Mizar occasionally appeared double and this way discovered the first double-lined spectroscopic 
binary star. Vogel measured radial velocities (RVs) of Algol and used them to prove that the known variations in the 
brightness of Algol are indeed caused by a "dark satellite revolving about it" (fVogcl 1890b). 

Vogel and collaborators built a series of prism spectrographs for the 30-cm refractor of the Potsdam Astrophysical 
Observatory IVogcl 1900). In 1888, they initiated a photographic radial- velocity program. Vogel was not the first 
to take a ph otograph of a stellar spectrum but improved the technique and obtained an RV precision of 2-4 km s~^ 
(|Vogellll891[) . Vogel (1906) and Pickering (1908) were both awarded the Bruce medal. Yet another Bruce medalist 
(1915), William W. Campbell (1862-1938) and collaborat ors carried out a larg e spectroscopic program at the Lick 
Observatory and discovered many spectroscopic binaries (jCampbell et al.lll911f ). They prepared the first catalog of 
spectroscopic binaries (Campbell & Curtis 1905). Vogel was the director of the Potsdam Observatory for 25 years, 
Pickering of the Harvard College Observatory for 42 years and Campbell of the Lick Observatory for 30 years. The 
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administrative duties must have been not too distracting those days. 

Both Vogel and Campbell recognized the importance of flexure and temperature control of a spectrograph (jVogell 
Il890d: 'Campbelll llSQSl . Il900f ) and Campbell also noted the issue of a slit illumination and its impact on RV precision 
ICampbcU 191^. Today a high RV precision is achieved by either using a very stable fiber fed spectrograph that 
is contained in a controlled environment (a vacuum tank) or using an absorption cell to superimpose a reference 
spectrum onto a stellar spectrum and this way measure and account for the systematic RV errors. Current state of 
the art precision is at the level of ~1 m s^^. It is however important to note that such a precision refers to single 
stars or at best single-lined spectroscopic binaries where the influence of the secondary spectrum can be neglected. In 
such a case, given a stable spectrograph, an RV measurement is essentially a measurement of a shift of an otherwise 
constant shape (spectrum). 

Radial velocities (RVs) of double-lined spectroscopic binary stars (SB2) can be used effectively to derive basic 
parameters of stars if the stars happen to be eclipsing or their astrometric relative orbit can be determined. It is quite 
surprising that the RV precision of double- lined binary stars on the ayerage has no t improved much over the last 100 
years (see Fig. [Ij. With the exception of our previous work (|Konackill20d5l[2009al) . the RV precision for such targets 
typically varies from ~0.1 km s~^ to ~1 km and clearly is much worse than what has been achieved for stars with 
planets or single-lined binary stars. The main problem with double-lined binary stars is that one has to deal with two 
sets of superimposed spectral lines whose corresponding radial velocities change considerably with typical amplitudes 
of 50-100 km s~^. In consequence a spectrum is highly variable and obviously one cannot measure RVs by noting a 
simple shift. 

We have developed a novel iodine cell based approach that employs a tomographic disentangling of the component 
spectra o f SB2s and allows one to mea sure RVs of the components of SB2s with a precision of the order 1-10 ms~^ 
([Kona cki 2009a; Konacki ct al.|[2"009bf ). Such quality RVs not only enable us to search for circumbinary extrasolar 
planets (Konacki ct al. 2009"bl) but also to determine basic parameters of stars with an unprecedented precision. In 
particular the masses of stars for non eclipsing SB2s can be easily determined with a fractional accuracy of the order 
of at least ^0.1% and often even ~0.01%. Moreover, we expect that the accuracy in masses will reach ~0.001% level 
when our method is applied to eclipsing binary stars. Such a level of precision is an order of magnitude higher than 
of the most accurate mass determination for a body outside the Solar System — the double neutron star system 
PSR B1913+16 (Nice et al. 2008). 

Below we present our precision RV data sets for five targets HD78418, HD123999, HD160922, HD200077 and 
IID210027 from our on-going TATOOINE (The Attempt To Observe Outer-planets In Non-single-stellar Environments) 
RV program to search for circumbinar y planets. All of them have been extensively observed with the Palomar Testbed 
Interferometer (jColavita et al.l[l999al PTI;). The archival PTI visibility measurements can be used to derive relative 
astrometric orbits of the binaries. These combined with our spectroscopic orbits allow for a complete orbital and 
physical description of the systems (with the exception of the radii of the components). In §2 we describe the RV 
measurements and their modeling and in §3 the visibility measurements and their modeling. In §4 we present the 
spectroscopic and astrometric orbital solutions and the resulting orbital and physical parameters of the binaries. A 
discussion is provided in §5. 

2. RADIAL VELOCITIES 

2.1. Iodine absorption cell and spectroscopic binary stars 

In the i odine cell (I2) techniq ue, the Doppler shift of a star spectrum AA is determined by solving the following 
equation (|Marcv fc Butledll992l i. 

I{\) = [F{\ + AX,)T{X + AXi,)](^PSF (1) 

where AAs is the shift of the star spectrum, AA/^ is the shift of the iodine transmission function T, (g) represents a 
convolution and PSF a spectrograph's point spread function (PSF). The parameters AAg, AA/^ as well as parameters 
describing the PSF are determined by performing a least-squares fit to the observed (through the iodine cell) spectrum 
/. For this purpose, one also needs (1) a high SNR star spectrum taken without the cell, F (the intrinsic stellar 
spectrum), which serves as a template for all the spectra observed through the cell and (2) the I2 transmission function 
T obtained with a Fourier Transform Spectrometer such as the one at the Kitt Peak National Observatory. The 
Doppler shift of a star spectrum is then given by AA — AXs — AA/j . Such an iodine technique can only be applied to 
single stars. This is dictated by the need to supply an observed template spectrum of a star in Eq. [TJ In the case of 
SB2s, it cannot be accomplished since their spectra are always composite and time variable. 

We can measure precise radial velocities of both components of an SB2 with an I2 absorption cell by performing the 
following steps. First, contrary to the standard approach for single stars, we always take two subsequent exposures of 
a binary target — one with and the other without the I2 cell. This way we obtain an instantaneous template which 
is used to model only the adjacent exposure taken with the cell. Next, we perform the usual least-squares fit and 
obtain the parameters described in Eq.[T] Obviously, the derived Doppler shift, AAi (where i denotes the epoch of the 
observation), carries no meaning since each time a different template is used (besides it describes a Doppler "shift" of 
a composite spectrum different at each epoch). However, the parameters (in particular the wavelength solution and 
the parameters describing PSF) are accurately determined and can be used to extract the star spectrum, /*'*(A), for 
each epoch t, by inverting the equation [1] 



r^*(A) = [/*(A) (^-1 PSF*]/T{X) 



(2) 
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where denotes deconvolution and PSF*, symbolically, the set of parameters describing PSF at the epoch t. Such 
a star spectrum has an accurate wavelength solution, is free of the I2 lines and the influence of a varying PSF. In the 
final step, the velocities of both compo nents of a b inary target can be measured with the well known two-dimensional 
cross-correlation technique TODCOR (iZucker fc M azch 1994) using as templates the synthetic spectra derived with 
the ATLAS 9 and ATLAS 12 programs C Kuruczlll99a and matched to the observed template spectrum, F{X). 

2.2. Iodine cell data pipeline 

In our data pipeline the reduction process involves the following procedures. An observed stellar spectrum, J, is a 
convolution of an intrinsic stellar spectrum, F, with a spectrographs' PSF, ^p. Such a convolution in a discrete form 
can be written as (jValenti et al.lll995t lEndl et al.ll2000D 

i-\-r 

h^Y. ^J^*-^' i = (3) 

where N is the number of pixels in the analyzed spectrum and r is the range of PSF such as t/'i = for |Z| > r. It is 
beneficial to work with an oversampled version of the above equation 

^ _ ^ ; = 1,...,7V (4) 

where o is the oversampling factor. This equation can be rewritten into two other useful formulas 

r / o— 1 \ 

(5) 





and 

r+0-1 / o-j-l \ 

[Y. ^APo+o^, ^ = l,...,iV (6) 

j=-r \ k=-j J 

They constitute a set of equations which can be written in the form I = RX where I represents a vector of the 
observed spectrum = 1, . . . , N and X represents a vector of the PSF = — r, . . . , r or of the intrinsic stellar 
spectrum Fk,k = 1, oN -t- 2r. The matrix R can be easily inferred from the appropriate sums above. By solving the 
equation 5 one can derive the PSF given the observed and intrinsic stellar spectrum and by solving the equation 6 one 
can derive the intrinsic stellar spectrum given the PSF and the observed spectrum. The latter process is obviously a 
deconvolution. The remaining needed set of equations is one that allows us to clean a stellar spectrum observed through 
an iodine cell from the iodine lines as it is symbolically described in the equation 2. This can be easily achieved by 
noting that in such a case ^1- in the equation 6 has to be replaced by F^T^. where Tfe is a known transmission function 
of the iodine cell. In consequence, the set of equations I = RX is replaced by I = RTX where T is a diagonal matrix 
with the diagonal elements equal to Tk and then solved for X {Xk — Fk, k — 1, oN + 2r). In our data pipeline the 
equations 5 and 6 are solved with the maximum entropy method using a commercially available MEMSYS5 software. 
Our data pipeline provides state of the art 1 m s~^ or better RV precision for single stars (see Fig. [J). 

In practice, the data reduction process is carried out as follows. Each observing night, a calibration exposure of 
a rapidly rotating B star or a quartz lamp is taken through an iodine cell. Such a spectrum in principle contains 
only the iodine lines and is used to determine the first wavelength solution and the first estimate of the PSF. This is 
done by taking F = 1 in the equation 1 and involves solving the equation 5 for tpk with Fk = Tk. The transmission 
function T of the iodine cell as a function of wavelength is obviously known. An observing sequence for an SB2 involves 
taking a pair of exposures one with and one without the I2 cell — the instantaneous composite template. The PSF 
is deconvolved from the template to obtain an intrinsic template spectrum, Fk, using the equation 6. The template 
is assumed to have a wavelength solution from the calibration exposure. Such a template is then used to determine 
the PSF and the wavelength solution for the exposure of SB2s taken with the cell. This is done by using the equation 
1 and involves solving the equation 5 where Fk is replaced by FkTk- In the final step the new PSF is used in the 
equation 6 to obtain a deconvolved and free of the I2 lines SB2 spectrum. This final spectrum is ready to be used with 
TODCOR as in .Konackl (,2005.1 . 

2.3. Tomographic disentangling of composite spectra 

The disentangling of the spectra of SB2s is a well known problem to which several more or less restricted solutions 
exist and ar e described in a rich lite rature. We have decided to follow the tomographic approach to the disentangling 
problem bv lBagnuolo &: GiesI (| 19911 ) which we numerically formulate and solve within the framework of MEMSYS5. 
The idea of the tomographic disentangling is presented in Fig. |3l Once the real template spectra for the components 
of a binary are available, it is no longer necessary to use TODCOR as one can simply use real templates and the 
original iodine equation [1] by replacing F(A -I- AA) with ^^(A -|- AA^) -|- SsS^{X + AA^) where 5^, 5^ represent the two 
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templates, AA"'^, AA'^ their shifts and 63 the brightness ratio. StiU TODCOR remains an important intermediate step 
to get the first approximate RVs that are subsequently used in the disentangling scheme. 

It is worth noting that the composite spectra have to be continuum normalized before they are tomographically 
disentangled and that the tomographic disentangling does not provide the real brightness ratio but only reproduces 
the depth of the spectral lines with respect to the normalized continuum of the composite spectrum. Hence ds is not 
the actual brightness ratio. Successful tomographic disentangling can be carried out based on as little as ~10 spectra 
given that they sample different radial velocities of the components. 

2.4. RVs and their errors 

The spectra used to derive RVs were collected with four telescopes and three spectrographs the Keck I/Hires, 
Shane/CAT/Hamspec and TNG/Sarg over the years 2003-2008. All three instruments are echelle spectrographs 
which in our setup provide spectra with a resolution of 67 000, 86 000 and 60,000 respectively and are equipped 
with iodine cells. The observing was carried out as described above. The typical signal to noise ratios (SNRs) per 
collapsed pixel of the collected spectra are ~250 for the Keck I/Hires, ^75-150 for the TNG/Sarg and ~50-150 for the 
Shane/CAT/Hamspec. These are for the composite observed spectra of SB2s. Hence for example a brightness ratio of 
2.3 (HD78418) corresponds to an SNR of ~175 for the primary component and an SNR of ~75 for the secondary. The 
RVs of the primaries are thus typically more accurate than of the secondaries. The formal RV errors were computed 
from the scatter between echelle orders. As it turns out these errors are underestimated. In order to obtain a reduced 
equal to 1 for a spectroscopic orbital solution and thus conservative estimates of the errors of the least-squares 
best-fit parameters, we add an additional error in quadrature to the formal errors (see §4). The most likely cause of the 
underestimation of the formal errors, in addition to the RV jitter of the stars, are subtle imperfections in the template 
component spectra obtained with the tomographic disentangling and used as a reference for RV computation. Figure[4] 
demonstrates how subtle the problem is. The sum of tomographically disentangled component spectra provide for an 
overall better match to observed spectra then a single composite observed spectrum. Still in many cases we do not yet 
reach a photon limited RV precision for SB2s. This problem is being investigated and future releases of more accurate 
RVs from our program are very likely. The current RVs used in this paper are listed in Table 1. 

2.5. RV modeling 

In a binary's center of mass coordinate system where the xy-plane is in the plane of the sky and the z-axis is directed 
away from the observer, the RVs may be modeled with the following equation 

Vtotal =Vr+Vz+ AVtides + AVsR GR (7) 

where Vz is a standard RV variation due to a Keplerian motion, AVudes is an RV variation due to tidal distortion of 
the components, AVsr gr is a relativistic contribution to RVs and Vr is the radial velocity of the center of mass of a 
binary (i.e. the gamma velocity). 

2.5.1. Keplerian motion and the light-time effect 
The standard Keplerian part of the RV variations is as follows. 

V = iV^,Vy,Vz) = ^ (-aP smE + aQ^/l-e^ coss) (8) 

n = 27r/P, r = a (1 — e cos £') , P = 1 coscj + m sino;, Q = — 1 sincj + mcosw, 
1 = (cos il, sin n, 0), m = (— cos i s'mfl, cos i cos ft, sini). 
In particular, Vz can also be traditionally expressed as 



I , p \ \ an sini f 1 + e E 

Vz ^ K {cos{ f + uj) + e cos uj) , K = =, tan - = W tan—, (9) 

VI -e^ 2 Vl-e 2 

where / is the true anomaly and E = E(t) is the eccentric anomaly given by the Kepler equation E — esinE = 
27r(t — Tp)/P, P is the orbital period, a, e, w, fi, Tp are the standard Keplerian elements — the semi-major axis of a 
component's orbit, the eccentricity, the inclination, the longitude of pericenter, the longitude of ascending node and 
the time of pericenter. The spectroscopic orbits of the components can be fully described using the two respective RV 
amplitudes Ki,2 = ai.2 n sin i/Vl — e^, with the longitudes of pericenter satisfying the relation (jJ2 = ij^i + tt and the 
remaining orbital parameters being the same. Note that here ai 2 refers to orbits of the components with respect to 
the center of mass of the binary. 

At the level of precision 1-10 m s^^ it is important to include the light-time effect within the binary's orbit. We 
achieve this by solving the following implicit equation 



r(ii,2) = ii,2 + ^i,2(ii,2)/c. 



^1.2(^1.2) =r- sini sin (/ + u;i.2) = 7- (l - e^)'/' sin (/ ^1,2) , (10) 

cn (1 + ecos/) ^ ' 
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for ti^2, where T is the observed moment of an RV measurement referred to SSB (i.e. corrected for the light-time 
effect within the Solar System), ti^2 is the actual moment and 1,2 refer to the primary and secondary respectively. 
Subsequently, ti^2 is used as the proper argument in the earlier equations. However, it is approximately true that 

tl,2 WT-Zi,2(tl,2)/C (11) 

and since .^i.2(^i. 2)/c is small, it can be eas ily derived that to the first order the RV contribution from the light-time 
effect is (see also lZucker fc Alexandei][2007[ ) 

AVlte,i,2 = -Kf 2 sin^ (/ + uji,2) (1 + e cos /) (12) 
c 

and can be used to estimate the amplitudes of AVlte.i,2, Kf2/c, which are 2.9 and 3.2 m s^^ for HD78418, 15.1 and 

16.0 m s-i for HD123999, 4.4 and 6.7 m s^^ for HD160922, 2.9 and 4.6 m s^^ for HD200077 and 7.8 and 20.2 m s^^ 
for HD210027, for the primary and secondary respectively. 

2.5.2. ndal effects 

The contributions of the tidally distorted binary components to the observed RVs were first noted by 'Sterne' fl941l). 
Such contributions arise from the non-spherical shapes of the binary's components and their non-uniform surface 
br ightness due to predom ina ntly gravitational and limb d arkenings. The corresponding RV variations were analyzed 
bvlWilson fc SofiaT (119761) and lKopal l|1980f ). iKopall ()1980n made an attempt to derive an analytic approximation while 
Wilson fc S ofia (1976) used a more numerical approach employing the Wilson-Davinney code (jWilson fc DevinnevI 
1971t rWilsoniil979l) for modeling eclipsing binary light curves to compute the "tidal" RV co ntribution. Recently, this 
effect was also analyzed by Eaton (2008). We follow the approach by Wilson & Sofia (|1976l ) and [Eaton (2008) where 
the "tidal " contribution to the RVs is modeled with the Wilson-Davinney code. To this end we use the light curve 
modeling program Ic which is part of the Wilson-Daviney code assuming the parameters of our targets as in Table 
3. The resulting RV variations are shown in Fig.[5l It should be noted that even though these effects are quite small 
they can have a significant qualitative effect if ignored by e.g. mimicking a small orbital eccentricity in an otherwise 
circular orbit. 



2.5.3. Relativtstic effects 

The relativistic description o f the motion of a binary star from the point of view of its radial velocities was derived 
by Kopeikin & Ozcrnov (1999"). In ou r model we adopt the following relativistic correction based on the equations 
76-86 from Kopeikin & Ozcrnoy (Il999[ ) 

AVsRGR = ^{V^ + V3) - - (^Vi + Ue) -—S-Vg^-(^ + Uss - (S • Vaf 
Zc c \Z / c c \ 2 

- is • Vg - is • Vg (V, us + Vy Ai„) {t - to) + (V, fis + Vy Hc.) [t - <o) + 

-Vr V, + -Vr (V, ^ls + Vy A*„) {t -to) + - (V, v., + Vy Vy) + - {\v^ + Uc\ (13) 

were Vr,Vt are the radial and tangential velocities of the center of mass of the binary, Ve,Ue are the velocity 
magnitude and gravitational potential of the observer with respect to the geocenter, Vg, Vq are the velocity vector 
and its magnitude of the geocenter with respect to the Solar System barycenter (SSB), Uss is the gravitational 
potential at the geocenter due to all major bodies in the Solar System, fiajfJ-s is the proper motion of the center of 
mass of the binary expressed in radians per time unit, Vx, Vx are the components of the tangential velocity vector such 
as 14 = pcua/ Vy = pc^s/ 1^ where n is the parallax, pc is one parsec, and S is the unit vector towards the center of 
mass of the binary. Vx, Vy, Vz as before denote the coordinates of the orbital velocity vector of a given component, V 
is its magnitude and Uc is the gravitational potential at the position of a binary component due to the gravitational 
field of its companion. Note that S varies slowly due to the proper motion vector /x, S = S -\- fj. [t — t^) and the term 
(Vx /i(S Vy ^a) {t — to) is not really relativistic but is obviously due to the varying S. The necessary quantities such as 
Vg, Ue, Uss, Vt,S we calculated using the JPL ephcmerides DE405 and the c atalogue posit i ons a nd proper motions 
of the targets with the help of the NOVAS library of astrometric subroutines bv lKaplan et al.l (|1989( ). 

The last term in the above equation is the combined effect of the transverse Doppler effect and the gravitational 
redshift and also the dominant term in the relativistic correction. It can be easily shown that its periodic part has the 
following form 



1,,2 \ _ r G'M2,i (Mi,2 2M: 



^v' + uc] =71,2 cos/, 71,2 = "'•'"^;;7" 'J7'^\. ^ 2, (14) 

2 JpertoMc Ca{Mi+M2) (1-6^) 

for the primary and secondary respectively where Mi, M2 are the masses of the components. One can employ the above 
equation to test GR by using 71^2 as a free parameter and fitting for it. However as noted by iKopeikin fc Ozernovl 
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(|1999f ). it is difficult in practice due to the coupling of this effect with the Keplerian part of the model. The above 
equation can be rewritten in the following form 

71,2 = ■ —. — e (15) 

c sm I 

and, as has been recently noted bv lZucker fc Alexander! (|2007T ). can in principle be used to derive the orbital inclination 
using just the RV measurements. 

Three of our targets HD78418, HD123999 and HD200077 have significant eccentricities and the resulting 71^2 are as 
follows: 4.8 and 5.9 m s~^, 9.7 and 10.1 m s~^, 8.1 and 11.0 m s~^ for the primary and secondary respectively. Only 
for HD78418, the procedure for deriving sin^ i proposed by IZucker fc Alexander! (|2007f) provided a value of 0.3832 
[i = 141.8 deg) close to the real one of 0.2985 {i = 146.9) as determined with the help of PTI astrometric data. 
However the error of sin^ i = 0.3832 is 1.8 so this is more of a coincidence. The method may yet turn out to be useful 
given sufficiently accurate and/or numerous RV data sets. 

Finally, let us note that our attempts to fit for lu which could be either due to the relativistic or tidal precession 
(jMazch 2008) have not produced any meanigful results. In four cases uj was at the level of its formal error. In one case 
(HD78418) UJ was technically relatively significant (2.4±0.6x 10~^ deg/day) but its inclusion in the model have not 
improved the best-fit rms. Hence, the orbital part of the description of the motion remains in its classic Newtonian 
form as in Eq. 8 throughout this paper. 

3. VISIBILITIES AND THEIR MODELING 

Interferometers such as PTI typically measure a fringe contrast - a normalized (by the total power received from 
a source) amplitude of the coherence function. The norma lized visibility of a binary star can be modeled with the 
following equation fsee lBodenll2000bl : iKonacki fc Landl2004[ ) 



^2 _ + r^Vi + 2rViV2 cos (27rB^ • As/A) 

(1 + 

where Vi , V2 are the visibilities of the components approximated with the visibility of a uniform disk of diameter 



^binary , \2 ^ ' 



2Ji{7r0B^X) 



2 



where B± = ||Bj^|| is the length of the projected baseline vector of a two aperture interferometer, r is the brightness 
ratio at the observing wavelength X {r — P2 /Pi where Pi , P2 are the total powers of the binary components at the 
given wavelength) and As = (A(5, Aa) is the separation vector between the primary and the secondary in the plane 

tangent to the sky. 

Th e separation vector between the primary and the secondary is given by the following equations (see e.g. IKovalevskvl 
[T995t iKamp 1967.1 

As(<) = {AS{t),Aa{t)) = Ka (p (cos E{t) ~e) + Q\/l- e? sin £;(t)) , (18) 

here a is the semi- major axis of the relative orbit and the first two coordinates of the vectors P, Q are used. Note 
that traditionally the North direction is the x-axis and the East direction is the y-axis for a coordinate system used in 
astrometry to model the orbital motion. In our figures showing the relative orbits, the x-axis is for the right ascension 
and the y-axis for the declination but the orbital motion is modeled as above. Note also that a relative astrometric 
orbit allows for a few possible solutions for the pair of angles namely cj, $7 as well as a; ± tt, f2 ± tt. The inclusion 
of a spectroscopic orbit allows one to determine the real w and traditionally is chosen to be less than tt. 

PTI provides a square of the normalized visibility amplitude, 1^^, as one can obtain its unbiased estimate (jColavital 
[l999b; Mozurkewich ct al. 1991). A normalized visibility of a point source should be by definition equal to 1. Since 
a real instrument is not perfect and does not operate in a perfect environment, this is typically not the case and the 
observed visibility is underestimated. In consequence each visibility measurement has to be calibrated. This is carried 
out by observing at least one calibrator in between the observations of a target. The calibrator is typically a single star 
which diameter is known and its vi sibility can be very well app roximated with the equation 17. Th en the calibrated 
visibility of a target is given by (Boden. Colavita. van Belle, fc Shao 1991 iMozurkewich et al.l[T99lh 

^calibrated ^rneasured /Kls (19) 

where 

^sys ^cal — measured/^cal — expected (^^) 

and V^a/_meQSMred ^ mcasurcd visibility and K-iz-ez^pected expected visibility of a calibrator given by the equation 
(9). 

Visibilities of our targets were extracted from the NASA Exoplanet Science Institute's (NExSci) database of the PTI 
measurements. They were subsequently calibrated using the excellent tools getCal (ver. 2.10.7) and wbCalib (ver. 
1.7.4) provided by NExSci. The data reduction was carried out using the default parameters of wbCalib. The software 
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provides {ti,V^^,ai,Xi,Ui,Vi} where U is the time of observation, V^^ is the calibrated visibility amplitude squared, 
(Ti its error, Ai is the mean wavelength for the observation, and {ui,Vi) are the components of the projected baseline 
vector. We do not list them as they can be easily obtained using the archival PTI data and the available tools. 

4. LEAST-SQUARES FITTING TO THE COMBINED RV/V^ DATA SETS 

We determined the best fit orbital parameters of the binaries in the standard way by minimizing the following 
function: 

N M M 

i—l i—1 i—1 

were V^^, (Ty2 , RVli, cTBVii a-nd RV2i, <jRV2i are respectively visibilities and their errors, radial velocities of the primary 
and their errors and radial velocities of the secondary and their errors; N and M denote the number of measurements. 
The symbols with a hat denote the respective model values computed as described in the previous sections. The fit was 
carried out using the Levenberg-Marquardt algorithm for solving the minimization problem which was incorporated 
into our own software for modeling and least-squares fitting of the visibilities and radial velocities. As mentioned 
above, the RV variations due to the tidal distortions of the components were modeled with the Wilson-Davinney's Ic 
code (version from 2007) which was also incorporated into our software. 

The least-squares fitting formalism allows one to compute the formal uncertainties (errors) of the best fit parameters. 
However such errors do not neccesarily correspond to the true uncertainties of the parameters. In order to provide 
conservative estimates of these errors we also accounted for the systematic errors in the modeling which are due to (1) 
the uncertainty in the projected baseline (w, v) (2) the uncertainty in the mean wavelength A (3) the uncertainties in 
the diameters, dj, of the calibrators (which affect the modeled visibility through the equation 20) and the diameters of 
the binary components (which affect the modeled visibility through the equation 17) and (4) the uncertainties in the 
parameters used to model the tidal RV effects. We assumed the following estimates for these uncertainties (1) 0.01 
percent in (u,w), (2) 0.5 percent in A and (3) 10 percent in the cahbrator and binary components diameters and (4) 
10 percent in all the parameters from Table 2 except for the temperatures for which we assumed an uncertainty of 1 
percent for HD210027 and 2 percent for the other targets; and for the mctallicities we assumed an uncertainty of 0.05 
dex. The errors of the parameters shown in Tables 4-6 include the contribution from the systematic errors. 

Our procedure of deriving RVs results in slightly different zero-points of the velocities for each data subset. These 
small differences are among others due to a mismatch, different for each component of a binary, of the synthetic 
templates and the real spectra. The synthetic templates, as already explained, are used in the first step of RV 
computations and the shifts are then carried over to the real disentangled component spectra. The other sources of 
the small zero-point differences are due to the use of different spectrographs and the accuracies with which the zero 
points of the template iodine cell transmission functions are known. The velocity shifts are determined together with 
the other parameters in the least-squares fits and are shown in Table 5. Even though the total formal errors of the 
gamma velocities are small, due to the above reasons it is unlikely that their realistic uncertainties are smaller than 
~0.1 km s-^ 

The full statistical detailes of the best-fit solutions are given in Table 5. For the RV data sets we determined the 
additional RV errors added in the quadrature to the formal RV errors by performing separate orbital fits to the RVs 
alone and finding such additional RV errors which provide ~ 1. For V'^ measurements we used the errors as 
computed with wbCalib. However by performing a preliminary fit to all the available V'^s and RVs, we found those 
V'^ measurements that significantly deviate from the best-fit solutions. Such deviant measurements are most likley 
due to poor weather conditions and/or an errant behaviour of the instrument. In the final fit we did not use them. In 
Figures 7,9,11,13,15 they are denoted with open (red) circles. The resulting total values are close or very close to 
1 which ensures one that the estimates of the errors of the parameters are realistic. 

4.1. Notes on the individual solutions 

The orbital parameters of the combined spectroscopic-astrometric solutions are shown in Tables 4-5 and the resulting 
physical parameters are shown in Table 6. The masses expressed in the solar masses were derived with GMq = 
1.3271244017987 x 10^" m^s^, a value used in the DE405 JPL ephemerides. Let us note that the official value of 
G recommended by CODATA (The Committee on Data for Science and Technology) is 6.67428 ± 0.00067 x 10^" 
jn3]jg-ig-2_ jjej^j^e Q jg known with a fractional error of 0.01% and in consequence the mass of the Sun and the masses 
of stars expressed in absolute units are known with only such precision. 

We also calculated limits to circuinbinary planets for three new systems HD123999, HD160922 and HD200077. They 
were calculated as in lKonacki et al.l ()2009br ) and are shown in Fig. 18. 

HD 78418 (75 Cnc, HE 3626, HIP 44892; V = 5.98 mag, K = 4-37 mag) is a ^19.4 d ay period binary with a spectral 
type G5IV-V. The most recent spectroscopic solution is bv lde Medeiros fc Udrvl ()1999f ) and is characterized by an RV 
rms of 461 m s~^. Our combined RV solution has an rms of 13.1 and 2 4.3 m for the primary and secondary 
respectively. The binary was resolved with PTI by iLane fc BodenI (|1999D but a detailed spectroscopic-astrometric 
analysis was never published. Using the PTI data archive, we extracted 440 V^ measurements spanning 1999-2008 
of which 415 were used in the final fit by adopting a cutoff at 0.1 for V^'s 0-C. The masses of the components are 
1.173±0.024 Mq (2.0% accuracy) and 1.011±0.021 Mq (2.1% accuracy) for the primary and secondary respectively. 
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This is a far more accurate mass determination than the one from lLane fc Boded ([19991). Th e masses and the absolute 
magnitudes in the K and H bands together with the Padova isochrones ( Marigo et al.l 12008?) allow us to estimate the 



age o f the system as about 2.5-4.0 Gyr. For the isochrones we a dopted a metallicity o f -0.09 (|Nordstrom "erall[2008l . 
\2004 ). The JHK photometry comes from the 2MASS catalogue ()Skrutskie et al.ll2006f l. 

HD 123999 (12 Boo, HR 5304, HIP 69226; V = 4.82 mag, K = 3.64) is a ^9 .6 day period binary with spectral 
type F8IV. The latest spectroscopic orbit was published bv lTomkin fc Fekell (|2006i) . It has an RV rms of 180 and 100 
m s^^ for the primary and secondary respectively. Our combined RV data sets are characterized by an rms of 34.0 and 
38. 3 m s~^ which is du e to somewhat wide s pectral lines of the s tars. The binary was resolved and studied with PTI 
bv lBoden et al.l ()2000af ) and most recently bv lBoden et al.l (120051 ) who obtained a combined spectroscopic- astrometric 
solution using the PTI measurements collected in the years 1998-2004. We extracted 374 measurements spanning 
1999-2008 and used 346 in the final fit (a cutoff at 0.09 for V^'s 0-C). The masses of the components are 1.4109±0.0028 
Mq (0.2% accuracy) and 1.3 677ib0.0028 M|7i (0.2 % accuracy) for the primary and secondary respectively. These values 
are co nsistent with those by iBoden et al.l ()2005D but about two times more accurate. The RVs used by iBoden et al.l 
(|2005[ ) were characterized by an rms of 400 and 490 m s~^. Assuming a solar metallicity as IBoden et al.l (j2005l ). the 
estimated age for the system is 2.5-2.9 Gyr. Note that as in IBoden et al.l (|2005| ). the secondary appears to be younger 
than the primary. Our smaller errorbars for the masses and absolute K and H band magnitudes make this discrepancy 
even more apparent. The JHK photometry comes from IBoden eFalj (|2000aD . 

HD 160922 (uj Dra, HR 6596, HIP 86201; V = 4.90 mag, K = 3.62 ma g ) is a ^5.3 day pe riod binary with a spectral 
type F5V. An improved spectroscopic solution was recently pubhshed bv lFekel et al.l (|2009( ). It is characterized by an 
RV rms of 190 m s~^. Our combined RV solution has an rms of 37.9 and 31.9 m s~^ for the primary and secondary 
respectively. This somewhat lower RV precision is due to relatively wide spectral lines of the components. The PTI 
archive provides 62 never published V^ measurements of which 61 were used in the final fit (V^'s 0-C cutoff at 0.08). 
We have obtained a spectroscopic-astrometric orbital solut i on wi th a small but statistically significant eccentricity 
(0.00220±0. 00031) which is in agreement with iFekel et al.l (|2009l ). However they have decided to adopt a circular 
solution afterall. The masses of the components are 1.46±0.16 Mq (11% accuracy) and 1.18±0.13 Mq (11% accuracy) 
for the primary and secondary respectively. The masses are derived for the first time for this system but their accuracy 
is poor due to a limited number of V^ measurements and the near face-on orbital configuration. Based on the derived 
masses and the K band absolute magnitudes we estimate the age of the system as 0.004-2.5 Gyr; a wide range due 
to large error bars both in th e masses and the K-band magnitudes. For the isochrones we adopted a solar metallicity 
(jNordstrom et al.l[2008L[200l . The JHK photometry comes from the 2MASS catalogue. 

HD 200077 (HIP 103641; V = 6.57 mag, K = 5.12 mag) is a ~112.5 day period bin ary with a spect r al typ e GOV. 
The spectrscopic orbital solution for both components was for the first time derived bv lGoldberg et al.l (|2002h . Their 
RVs are characterized by an rms of 600 and 2150 m s"^ for the primary and secondary respectively. Our combined 
RV solution has an rms of respectively 8.5 and 30.4 m s^^. The binary was resolved with PTI but this has never 
been published before. The PTI archive provides 395 V^ measurements spanning 2001-2007 of which 329 were used 
in the final fit (V^'s 0-C cutoff at 0.09). The masses of the components derived for the first time are 1.1860±0.0057 
Mq (0.48% accuracy) and 0.9407±0.0049 Mq (0.52% accuracy) for the primary and secondary respectively. Based on 
the masses, K and H band absolute magnitudes the estimated age is 2.9-4.3 Gyr. For the isochrones we adopted a 
metallicity of -0.14 (jNordstrom et al. 2008, 2004). The JHK photometry comes from the 2MASS catalogue. 

HD 210027 (l Peg, 24 Peg, HR 8430, HIP 109176; V = 3.76 mag, K = 2.56) is a ~10.2 day period binary with 
spectral types F5V/G8V. It is one of the first SB2s resolved with PTI. The most recent orbital solution is a combined 
spectroscopic-astrometric soluti on by IBoden et al.l (1999"). It is based on the PTI V^ measurements from 1997 and the 
archival RV measurements from iFekel fc TomkinI (l983) which are characterized by an rms of about 600 and 700 m s 
for the primary and secondary. Our combined RV solution has an rms of respectively 17.1 and 69.7 m s^^. The PTI 
archive provides 299 V^ measurements spanning 1997-2006 of which 266 were used in the final fit (V^'s 0-C cutoff at 
0.088). The masses of the components are 1.33249±0.00086 Mq (0.065% accuracy) and 0.83050±0.00055 Mq (0.066% 
accuracy) for the prim ary and secondary respectively. This is about a 20 times more accurate mass determination 
than the one from .Boden et al.l ( (199 9) and constitues the most accurate mass measurement for a normal star. Based 
on the mas s es, K and H band absolute magnitudes the estimated age is 4-663 Myr w hich i sin agreement with 
iMorel et al.l (120001). For the iso chrones we adopted a solar metalHcity (jMorel et al.ll200Q) . The JHK photometry 
comes from iBouchet et all ()1991t ). Note that our orbital solution has a small but statistically significant eccentricity 
(0.001764±0. 000063). This combined with the systems' young age may prove useful! in studying the tidal history of 
HD 210027. 

5. DISCUSSION 

There are a few ways to determine accurate masses of normal stars. Perhaps the most classic one is through absolute 
astrometric orbits of both components of a binary. However this constitues a challenging measurement and in practice 
masses are typically derived using diverse data sets e.g. (1) by combining RVs, relative astrometry and parallax for 
single- lined spectroscopic binaries, (2) as in this paper by combining relative astrometry and RVs for double-lined 
spectroscopic binaries and (3) by combining light curves and RVs for eclipsing double-lined binaries. The last method 
is the most useful one as it not only provides the most accurate masses due to a convenient edge-on geometry (note 
that the masses are d erived from their res pective M sin^ i) but also enables one to determine the radii of stars. 

In a recent review, iTorres et al.l (|2009|) collect 118 detached binary stars (including 94 eclipsing) with the most 



High precision orbital and physical parameters of SB2s 



9 



accurate mass determinations in the literature. These are denoted with open circles in Fig. 18. Even though our 
targets are not eclipsing and the orbital inclinations and their errors have a significant impact on the precision in 
masses, for two s tars HD1239 99 and HD210027 we have obtained more accurate mass determinations than for any 
of the stars from iTorres et al.l (j2009) . The accuracies of the masses of these two binaries are in the precision range 
covered by only close double neutron star systems characterized with radio pulsar timing. In particular, the masses 
for HD210027 riv al in precision th e mass determination of the components of the relativistic double pulsar system 
PSR J0737-3039 ([Nice et al.l 120081 ). If our targets were all eclipsing and the accuracy in masses was limited by our 
RVs alone, the accuracy would be in the range 0.02% to 0.42% (the fractional accuracy of M sin^i). The lower limit 
of this range is equal to the mas s accuracy of PS R B1913+16 which has the most accurate mass determination for a 
body outside the Solar System (jNice et al.ll2008f l. In fact if we adopt our RV precision and use the achievable from 
the ground precision in the orbital inclination angle for eclipsing binaries, we can expect to obtain masses with a 
fractional precision of ~0.001% (see Fig. 18). Clearly, our RV technique for double-lined and echpsing binary stars 
opens an exciting opportunity for derving masses of stars (and other parameters) with an unprecedented precision. 
These combined with parallax measurements from e.g. the planned GAIA astrometric mission and hopefully accurate 
abundance determinations should produce an outstanding set of parameters to tests models of the stellar structure 
and evolution. 

We thank the California and Carnegie Exoplanet Search team, and Geoff Marcy in particular, for allowing us access 
to their precision velocimetry tools at Lick Observatory. This work benefits from the efforts of the PTI collaboration 
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PTI's night assistant Kevin Rykoski for his efforts to maintain PTI in excellent condition and operating PTI. This 
research has made use of the Simbad database, operated at CDS, Strasbourg, France. MWM acknowledges support 
from the Townes Fellowship Program. M.K. is supported by the Foundation for Polish Science through a FOCUS grant 
and fellowship and by the Polish Ministry of Science and Higher Education through grant N203 005 32/0449. Part of 
the algorithms used in this analysis were developed during the SIM Double Blind Test, under JPL contract 1336910. 
This research has made use of the Simbad database, operated at CDS, Strasbourg, France. The observations on the 
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collaboration supported by the Research Infrastructures Programme of the European Commissions Sixth Framework 
Programme. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint 
project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of 
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Fig. 1. — Radial velocity precision (an rms from the orbital fit) for the primaries of double-line d spe ctroscopic binaries as a function of 
the publication date based on the ninth catalogue of spectroscopic binary orbits ( Pourbaix et al , 2 0041) . It is worth noting that already in 
the early XX century an RV precision of several km s~^ for double-lined binaries was possible l|Plummer et al.lll908l . HD83808). Current 
and previous precision range from our method is delimited with the red stars (note the logarithmic scale). 
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Fig. 2. — A series of the precision RV measurements of a single star HD185144 (V=4.7 mag, KOV) taken with the Keck I/HIRES and 
reduced using our iodine cell data pipeline. The data were taken on 5 nights: two in July (5 consecutive spectra each night) and three 
in October 2005 (four consecutive spectra each night). The average SNR was 600 for the template exposure as well as for all the spectra 
taken with the iodine cell. 
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Fig. 3. — An SB2 as nature's realization of the tomographic imaging. The composite spectrum of an SB2 taken at any orbital phase can 
be interpreted as a result of imaging of a two-layered object. By observing an SB2 at several different orbital phases and hence different RVs 
of its components, one can carry out an equivalent of tomographic imaging and eventually be able to disentangle the component ("layer") 
spectra. 
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Fig. 4. — Left and right (a): a piece of a spectrum of HD4676 as a function of wavelength as seen through an iodine cell. Left (b): a 
template composite spectrum of HD4676 taken just after the iodine cell exposure shown in (a). Right (b): the corresponding continuum 
normalized disentangled component spectra of HD4676 shifted accordingly to their respective RVs. Left and right (c): the residuals after 
subtracting (b) from (a) and accounting for the iodine cell lines. The residuals from the disentangled component spectra are smaller than 
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Fig. 11. — Visibility measurements of HD160922 as a function of time (a), their best-fit residuals as a function of time (b) and histogram 
(c). The measurements used to determine the best-fit orbital solution are denoted with filled circles. The corresponding orbital coverage 
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Fig. 12. — Observed and modelled radial velocities of HD200077 as a function of the orbital phase (a), their best-fit residuals as a function 
of the orbital phase (b) and time (c). The histograms of the residuals for the primary and secondary (d). The Keck I/Hires is denoted 
with circles, Shane/CAT/Hamspec with triangles and TNG/Sarg with stars. 
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Fig. 13. — Visibility measurements of HD200077 as a function of time (a), their best-fit residuals as a function of time (b) and histogram 
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Fig. 15. — Visibility measurements of HD210027 as a function of time (a), their best-fit residuals as a function of time (b) and histogram 
(c). The measurements used to determine the best-fit orbital solution are denoted with filled circles. The corresponding orbital coverage 
and the relative orbit is shown in the panel (d). 
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Fig. 16.— The isochrones in the mass - absolute magnitude plane for HD78418 (top), HD123999 (middle) and HD160922 (bottom) 
together with the corresponding best-fitting isochrones. 
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Fig. 17. — The isochrones in the mass - absolute magnitude plane for HD200077 (top two panels) and HD210027 (bottom four panels) 
together with the corresponding best-fitting isochrones. In the case of HD210027 their companions are shown separately so as it is easier 
to compare the errorbars with the isochrones. 
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Fig. 19. — Fractional accuracy in masses for the binary stars with the most accurate mass determinations in the hterature. Open circles 
denote the binaries from a recent review by Torres ct al. (2009). Filled (blue) circles denote the double neut ron star system s B1913+16, 
B2127+11C, B1534+12, J0737-3039, J1756-2251 and J1906+0746 characterized with radio pulsar timing (N ice et al.l 120081) . Filled (red) 
rectangles denote the masses for HD78418, HD123999, HD200077 and HD210027 from this paper and a mass determination for AI Phe 
from our recent work (Helminiak et al.ii2009, ). The hatched area is the expected precision in masses for double-lined eclipsing binary stars 
assuming a mass range of the components of 0.5-1.5 Mq, orbital period range of 3-23 days, orbital inclination range of 85-90°, inclination's 
error range of 0.05-0.3° and radial velocity amplitudes' error range of 1-31 m s~^. 
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TABLE 1 

Radial velocities of HD78418, HD123999, HD160922, HD200077 and HD210027 



Target 


Time 
(TDB- 2400000.5) 


RVi 

(km s~-^) 


(km s-i) 


O-Ci 
(km s-i) 


(km s-i) 


RVa 

(km s"-*-) 


0-2 

(km s-i) 


O-Ca 
(km s"-"-) 


£2 


Inst. 


HD78418 


53094.347119 
53094.395761 


-14.83066 
-14.91605 


0.00854 
0.01039 


0.00538 
0.00241 


0.00299 
0.00663 


38.51581 
38.60908 


0.01243 
0.01267 


0.00234 
-0.00023 


0.00580 
0.00628 


HO 
HO 




54101.570433 
53744.479316 


-7.42109 
34.57326 


0.00876 
0.00817 


-0.00731 
0.00352 


0.00357 
0.00167 


29.91678 
-18.81645 


0.01302 
0.01216 


0.01759 
0.01278 


0.00697 
0.00518 


HN 
HN 




54788.464774 
54430.494032 


18.30307 
-2.05150 


0.00696 
0.00974 


-0.00953 
-0.00216 


0.00485 
0.00836 


0.16661 
23.77923 


0.03029 
0.03049 


0.01501 
0.01071 


0.01374 
0.01416 


H 
H 




54190.964445 
54191.956955 


34.26932 
29.78440 


0.01290 
0.01249 


0.00397 
-0.01327 


0.01290 
0.01249 


-19.07279 
-13.94505 


0.01534 
0.01050 


0.03793 
-0.01623 


0.01534 
0.01050 


S 
S 


HD123999 


53454.573974 
53456.545316 


-51.41296 
39.94170 


0.03155 
0.03231 


0.00976 
-0.02204 


0.00585 
0.00911 


72.66300 
-21.65504 


0.03675 
0.03623 


0.00268 
-0.04321 


0.01119 
0.00937 


HN 

HN 




53978.858598 
54247.021370 


24.88006 
48.74097 


0.04693 
0.03697 


0.00628 
0.07561 


0.03523 
0.02015 


-6.43652 
-31.00316 


0.03679 
0.03535 


-0.01637 
-0.04846 


0.03679 
0.03535 


S 
S 




54107.560803 
54108.572117 


-52.70193 
-25.53836 


0.02118 
0.04058 


-0.02634 
-0.01830 


0.01388 
0.03730 


73.95953 
45.98819 


0.03339 
0.05134 


-0.02783 
0.02694 


0.01819 
0.04303 


H 
H 



Note. — 0-12 denote the total errors used in the least-squares fits and €1^2 denote the internal errors. The last column denotes 
the spectrograph used to obtain a measurement. HN stands for Keck I/Hires with the upgraded detector, HO with the old detector, S 
stands for TNG/Sarg and H for Shane/CAT/Hamspec. The rest of this data set and the other RV sets are available via the link to the 
machine-readable version. 



TABLE 2 

Assumed Parameters for HD78418, HD123999, HD160922, HD200077 and HD210027'' 



Parameter HD78418 HD123999 HD160922 HD200077 HD210027 



Eff. temperature, primary, Ti (K) 


6000 


6130 


6500 


6000 


6642 


Efl. temperature, secondary, T2 (K) 


5900 


6230 


5900 


5500 


4991 


Potential, primary, f2i 


27.1 


12.03 


14.3 


26.7 


20.0 


Potential, secondary, O2 


34.9 


15.21 


14.1 


26.5 


22.5 


Synchronization factor, primary, Fi 


1.51 


1.5 


1.0 


6.57 


1.0 


Synchronization factor, secondary, F2 


1.51 


1.5 


1.0 


6.57 


1.0 


Gravity darkening exponent, primary, gi . 


0.4 


0.35 


0.3 


0.35 


0.3 


Gravity darkening exponent, secondary, 92 


0.4 


0.35 


0.3 


0.4 


0.4 


Albedo, primary, Ai 


0.5 


0.5 


0.5 


0.5 


0.5 


Albedo, secondary, A2 


0.5 


0.5 


0.5 


0.5 


0.5 


Metallicity 


-0.09 


0.0 


0.0 


-0.14 


0.0 


Apparent diameter, primary, 61 (mas) .... 


0.45 


0.638 


0.5 


0.26 


1.06 


Apparent diameter, secondary, 82 (mas) . . 


0.30 


0.480 


0.4 


0.21 


0.6 



The more accurate parameters for HD123999 and HD210027 come from respectively jBoSeneFaD l l2005l) and IMorel et al.l (12000) . The 
remaining ones are our best estimates based on the available information about the targets. 



TABLE 3 

Calibration stars for visibility measurements of HD78418, HD123999, HD160922, HD200077 and HD210027^ 



Target 


Oalibrator 


Spectral 


Magnitude 


Angular Separation 


Apparent 






Type 




irom larget (uegj 


Diameter (mas) 


HD78418 














HD79452 


G6III 


6.0 V, 3.8 K 


8.1 


0.79 




HD73192 


K2III 


6.0 V, 3.3 K 


9.0 


1.38 


HD123999 














HD121107 


G5III 


.5.7 V, 3.6 K 


8.2 


0.70 




HD128167 


F3V 


4.5 V, 3.6 K 


7.1 


0.79 




HD123612 


K5III 


6.6 V, 3.0 K 


0.9 


1.29 


HD160922 














HD154633 


G5V 


6.1 V, 4.5 K 


5.4 


0.31 




HD158633 


KOV 


6.4 V, 4.4 K 


1.8 


0.59 




HD168151 


F5V 


5.0 V, 3.9 K 


5.7 


0.56 


HD200077 














HD 192640 


A2V 


4.9 V, 4.9 K 


9.5 


0.71 




HD192985 


F5V 


5.9 V, 4.8 K 


9.6 


0.39 




HD199763 


G9III 


6.5 V, 4.3 K 


9.9 


0.78 


HD210027 














HD211006 


K2III 


5.9 V, 3.2 K 


3.6 


1.06 




HD211432 


G9III 


6.4 V, 4.2 K 


3.2 


0.70 




HD215510 


G6III 


6.3 V, 4.2 K 


10.7 


0.85 




HD210459 


F5III 


4.3 V, 2.5 K 


7.9 


0.81 



The adopted diameters of the calibrators are determined from their effective temperature and bolometric flux derived from archival photometry using getCal tool (ver. 2.10.7) 
supplied by the NASA Exoplanet Science Institute. 



TABLE 4 

Best-fit Orbital Solutions for HD78418, HD123999, HD160922, HD200077 and HD210027'' 



Parameter 


HD78418 


HD123999 


HD160922 


HD200077 


HD210027 


Apparent semi-major axis, a (mas) 


5.8696(96) 


3.4706(55) 


3.469(17) 


14.453(18) 


10.329(16) 


Period, P (d) 


19.412347(23) 


9.6045601(36) 


5.2797766(44) 


112.5131(14) 


10.2130253(16) 


Time of periastron, Tp (TDB-2400000.5) 


53895.4025(24) 


54099.93572(70) 


54348.583(83) 


53830.168(14) 


52997.378(52) 


Eccentricity, e 


0.19494(11) 


0.19214(15) 


0.00220(31) 


0.66227(51) 


0.001764(63) 


Longitude of the periastron, lo (deg) 


283.389(39) 


286.832(29) 


314.8(5.6) 


197.071(25) 


272.8(1.8) 


Longitude of the ascending node, Cl (deg) 


171.892(85) 


80.49(10) 


1.23(32) 


89.403(28) 


176.262(75) 


Inclination, i (deg) 


146.88(25) 


107.95(12) 


151.4(1.1) 


118.681(80) 


95.83(12) 


Magnitude difference (K band), AK 


1.1445(131) 


0.601(13) 


0.841(18) 


1.1968(88) 


1.675(15) 


Magnitude difference (H band), AH 


1.1726(349) 


0.667(30) 




1.263(38) 


1.75(11) 


Velocity amplitude of the primary, Ki (km s"'^) . 


26.4961(35) 


67.189(11) 


36.254(16) 


29.373(82) 


48.4757(39) 


Velocity amplitude of the secondary, K2 (km s~^) 


30.7579(65) 


69.311(14) 


44.720(16) 


37.03(11) 


77.777(16) 


Gamma velocity, Vr (km s~^) 


9.7478(60) 


9.646(13) 


-14.011(19) 


-36.009(18) 


-4.6504(34) 


^ The numbers in parentheses are the 1-sigma errors 


in the last digits quoted. 











TABLE 5 

Various Fit Parameters for HD78418, HD123999, HD160922, HD200077 and HD210027'' 



Parameter 

Velocity offsets 

Secondary vs primary (km s""*^) 

Hires old vs new detector, primary, vii (km s~^) 

Hires old vs new detector, secondary, vi2 (km 

Hamspec vs Hires new detector, primary, V21 (km s~^) 

Haiuspcc vs Hires new detector, secondary, V22 (km s~^) 

Sarg vs Hires new detector, primary, V31 (km s~^) 

Sarg vs Hires new detector, secondary, U32 (km s"'^) 

Hamspec vs Sarg, primary, V41 (km s~^) 

Hamspec vs Sarg, secondary, V42 (km s~^) 

Least-squares fit parameters 

Number of RV measurements, total 

Number of RV measurements, Keck/Hires 

Number of RV measurements, Shane/CAT/Hamspec 

Number of RV measurements, TNG/Sarg 

Number of V^ measurements 

Additional RV error, Keck/Hires, primary /secondary (m s^-'^) 

Additional RV error, Shane/CAT/Hamspec, primary/secondary (m s~^) 

Additional RV error, TNG/Sarg, primary /secondary (m s~^) 

Combined RV rms, primary/secondary (m s^^) 

Keck/Hires RV rms, primary /secondary (m s^-*^) 

Shane/CAT/Hamspec RV rms, primary/secondary (m s~^) 

TNG/Sarg RV rms, primary /secondary (m s~^) 

V^ rms 

RV 1 primary/secondary 

y^x^ 

Degrees of freedom, DOF 

Total reduced x^, X^/DOF 



HD78418 


HD123999 


HD160922 


HD200077 


HD210027 


0.239(12) 


0.032(21) 


0.116(29) 


0.340(42) 


0.333(16) 


-0.0009(71) 








-0.0053(38) 


-0.0131(11) 








0.003(21) 


0.0358(67) 


-0.010(15) 






-0.0141(68) 


0.054(12) 


0.042(18) 






0.026(27) 


-0 3059(87~1 


-0 182('24') 




-0 3092('79~l 


-0 3026l'93') 


-0.286(10) 


-0.151(25) 




-0.229(23) 


-0.256(32) 






0.262(21) 










0.344(21) 






58 


64 


36 


26 


144 


26 


16 




18 


102 


24 


34 


20 




30 


8 


14 


16 


8 


12 


415 


346 


61 


329 


266 


8.0/11.0 


31.0/35.0 




5.0/26.0 


10.5/65.0 


5.0/27.0 


16.0/28.0 


27.0/15.0 




18.5/68.0 


0.0/0.0 


31.0/0.0 


37.0/25.0 


0.0/26.0 


10.5/38.0 


13.1/24.3 


34.0/.38.3 


37.9/31.9 


8.5/30.4 


17.1/69.7 


6.9/14.3 


30.2/35.8 




6.7/17.1 


14.1/64.6 


18.2/33.8 


36.3/37.6 


25.6/14.8 




22.7/84.2 


8.2/17.4 


40.2/26.8 


35.9/37.1 


10.4/25.9 


12.3/55.8 


0.0351 


0.0261 


0.0274 


0.0381 


0.0281 


31.60/32.99 


28.21/28.20 


19.31/18.19 


11.28/13.41 


70.40/64.94 


433.71 


482.46 


50.59 


522.01 


248.29 


454 


393 


83 


355 


391 


1.096 


1.371 


1.061 


1.608 


0.981 



^ The numbers in parentheses are the 1-sigma errors in the last digits quoted. 
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Physical Parameters for HD78418, 



TABLE 6 
HD123999, HD160922, 



HD200077 AND HD210027'' 



Star 



Pcirameter 

Semi-major axis, 01^2 (AU) 

M sin^i (M0) .' 

Mass, M (Mq) 

Mk,2MASS (mag) 

Mh,2MASS (mag) 

Parallax, k (mas) 

Distance, d (pc) 

Semi-major axis, ai_2 (AU) 

M sin^ j (Mq) 

Mass, M (Mq) 

Mk,2MASs (mag) 

^H,2MASS (mag) 

Parallax, k (mas) 

Distance, d (pc) 

Semi-major axis, ai_2 (AU) 

M sin^j (Mq) 

Mass, M (Mq) 

Mk,2MASS (mag) 

Mh,2MASS (mag) 

Parallax, k (mas) 

Distance, d (pc) 

Semi-major axis, ai_2 (AU) 

M sin^i (Mq) .' 

Mass, M (Mq) 

Mk,2MASS (mag) 

Mh,2MASS (mag) 

Parallax, k (mas) 

Distance, d (pc) 

Semi-major axis, 01^2 (AU) 

M sin^ j (Mq) .' 

Mass, M (Mq) 

Mk,2MASS (mag) 

Mh,2MASS (mag) 

Parallax, k (mas) 

Distance, d (pc) 



Primary 



Secondary 



HD78418 



HD123999 



HD 160922 



HD200077 



HD210027 



0.084874(11) 
0.191350(50) 
1.173(24) 
2.223(16) 
2.397(37) 



0.061193(10) 
1.21467(33) 
1.4109(28) 
1.254(18) 
1.314(20) 



0.036728(16) 
0.160409(88) 
1.46(16) 
2.16(25) 



0.25945(74) 

0.8008(34) 
1.1860(57) 
2.389(24) 
2.412(23) 



0.0457446(37) 
1.31190(29) 
1.33249(86) 
2.5125(70) 
2.606(19) 



0.098526(21) 
0.164836(45) 
1.011(21) 
3.367(18) 
3.570(45) 

32.004(52) 
31.246(51) 

0.063125(13) 
1.17749(33) 
1.3677(28) 
1.855(19) 
1.981(26) 

27.917(44) 
35.820(57) 

0.045306(17) 
0.130039(80) 
1.18(13) 
3.00(25) 

42.29(27) 
23.65(15) 

0.32711(96) 

0.6352(30) 
0.9407(49) 
3.585(25) 
3.675(36) 

24.640(32) 
40.571(52) 

0.073395(15) 
0.81767(22) 
0.83050(55) 
4.187(14) 
4.353(89) 

86.70(14) 
11.534(19) 



^ The numbers in parentheses are the 1-sigma errors in the last digits quoted. 



